Extended staggered-flux phases in two-dimensional lattices 
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Based on the so called t-(j> model in two-dimensional (2D) lattices, we investigate the stabilities of 
a class of extended staggered-flux (SF) phases (which are the extensions of the y/2 x \/2 SF phase 
to generalized spatial periods) against the Fermi-liquid phase. Surprisingly, when away from the 
nesting electron filling, some extended-SF phases take over the dominant SF phase (the v2 x \[2 
SF phase for the square lattice, a 1 x \/3 SF phase for the triangular one), compete with the 
Fermi-liquid phase in nontrivial patterns, and still occupy significant space in the phase diagram 
through the advantage in the total electronic kinetic energies. The results can be termed as the 
generalized Perierls orbital-antiferromagnetic instabilities of the Fermi-liquid phase in 2D lattice- 
electron models. 

PACS numbers: 71.70.-d, 75.10.-b, 75.10.Lp, 71.10.Hf 



I. INTRODUCTION 

An intriguing orbital-current-carrying state in the 
square lattice, the staggered-flux (SF) phase 1 (also 
known as orbital antiferromagne^, d-density wave 3 state, 
or charge flux phase) which has the spatial period \[2 x 
has been the focus of consideration for a long time. 
This y/2 x \/2 SF phase was first considered in excitonic 
insulators^, then in high-T c superconductors 1 *. The ab- 
sence of experimental vindications made it discarded in 
favor of other conventional orders. Recently, it regains 
new attentions since unusual experimental findings in 
two systems: a hidden order in the heavy-fermion com- 
pound URu2Si2^; a pseudogap phenomena (See the re- 
view by Timusk et aL 5 ) in the underdoped region of high- 
Tc cuprates 3 ". 

In parallel with the experimental contexts, numeri- 
cal signatures in the t-J model&i and evidences in two- 
leg ladder models (hall-filled t-U-V-M, doped t-J-V- 
V and t-J±-U-V±—) of its existence have been found. 
It is also being discussed now in many new theoreti- 
cal contexts: the double-exchange modelifi of colossal- 
magnetoresistance manganites; the repulsive SU(iV) 
Hubbard Modelii of ultracold fermionic atoms in opti- 
cal lattices; models with ring-exchange interactions 1 ^. 

For the triangular lattice 1 x \/3 SF phase with 

7r/2-nux per plaquette has been considered as the mean- 
field ansatai 3 - for the Heisenberg model. Recently found 
layered sodium cobalt oxide system NazCoO^j/H^O 14 
makes it a reasonable reference state for a range of elec- 
tron fillingiii. 

To the best of our knowledge, the formal exploration 
of other possible staggered-flux phases has never been 
made in previous studies (except a suggestion in study- 
ing the double-exchange modeliS). We would address the 
following problems: can a class of extended-SF phases 
with generalized spatial periods, which are still macro- 
scopically space- inversion and time-reversal invariant, be 



stabilized in a 2D lattice-electron model? if so, how do 
these extended-SF phases compete with each other and 
with the Fermi-liquid phase? is this kind of stabilization 
a generic feature of 2D lattice-electron models? 

In this report, we try to answer the above questions 
through numerical studies of a class of extended-SF 
phases in the so-called t-<f> modelii which is the sim- 
plest one being able to generate the \/2 x V2 SF phase 
spontaneously. This model has also been shownil to be 
closely related to the large-A^ SU(iV) t-J model (where 
the y/2 x y/2 SF phase is first realized as a ground stated) 
and the weak-coupling Hubbard model. If the idea about 
extended-SF phases can be realized in this simple model, 
then it should be further taken into account and be tested 
in the future studies of other more realistic models. 

For each extended-SF phase, we first investigate the 
variation of the total kinetic energy (TKE) of tight- 
binding electrons (electron-hopping part of the t-<p 
model) versus the flux parameter <f>. Further, includ- 
ing the "magnetic" energy of the flux itself, quantum 
phase transitions among these extended-SF phases and 
the Fermi-liquid phase are discussed in details through 
quantum phase diagrams. 



II. NOTATION AND FORMULATION 

The notation for an extended-SF phase is generally 
written as SF^ XS , where L=S,T represent the square lat- 
tice and the triangular one, respectively, and r x s rep- 
resents the spatial period (i.e., the size of a selected unit 
cell). Sometimes a Greek letter is added in the notation 
to distinguish among the extended-SF phases which have 
the same spatial period(e.<?. SF^ x ^ q and SF^^ ). 
We will consider 8 kinds of extended-SF phases in each 
lattice. The flux configurations of various extended- 
SF phases are illustrated in Fig. ^ These extended-SF 
phases are chosen by examining their symmetries and 
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simplicities. More complicated extended-SF phases can 
also be obtained by the generalized routines, however, 
they have less importance both theoretically and experi- 
mentally. 
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FIG. 1: Illustration of extended-SF phases in 2D lat- 
tices. Shaded plaquette has a flux <f>, and the blank 
one has a flux — < 
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The following U(l) gauge-invariant tight-binding 
Hamiltonian describes the t-6 modeli^: 



H 



£(< 



He 



p 



(1) 



Here (/>ij is in units of 0o/27r (0o = hc/e is the flux quan- 
tum). The nearest-neighbor hopping integral of electrons 
t is modified as tij = iexp(z^y ) due to the Aharonov- 
Bohm effect. The second term of Eq. gives the "mag- 
netic" energy. <j) v is the sum of 0y (clockwisely) along 
a plaquette p. To capture the major physics, the spatial 
fluctuation of \4> p \ is neglected. 

For each extended-SF phase, after taking a specific 
gauge (physical quantities are gauge-invariant as demon- 
strated beforei^), the first part of Eq. Q can be con- 
verted to fc-space and then be directly diagonalized, leads 
to several subbands of the energy spectrum. The sym- 
metries of the extended-SF phases in Fig. ^ ensure that 
the number of subbands are not too large. Among these 
extended-SF phases, the energy spectra of three simpler 
ones, SF^ J&1&, SFj xl and SF^ M i&, have been 

' V2x V2 lxV3,a ' 

clarified in details before. 



After getting the energy spectrum (numerically for 
large Hamiltonian matrices), the total density of states 
(DOS) D(uj) can be obtained by summing up the contri- 
butions from all subbands. Thus the TKE per site for 
electron-filling factor v at zero temperature is given by 



b,v) 



cLjD(uj)u>. 



(2) 



And the chemical potential fi is determined through 

dwD(uj). (3) 



Where oj m i n is the lower limit of the energy spectrum. 
Both the DOS and TKE are invariant under the trans- 
formations <fi — > —<f> and <p — > 2ir — <fi. Therefore, it is 
enough to consider Ekin for < 4> < it. For = (the 
Fermi-liquid phase) or ir (also called the 7r-flux phased in 
the square lattice), all extended-SF phases in a 2D lattice 
are equivalent and hence have the same Ekin ■ 

Then through the self-consistent determination of the 
flux parameter 4> by minimizing the total energy per site 
determined by the t-cf> model (TKE plus the "magnetic" 
energy), we can plot the phase diagram in the v-K pa- 
rameter space by labelling the lowest-energy phase (an 
extended-SF phase or the Fermi-liquid phase) in the v- 
K parameter space. Therefore for a given electron-filling 
factor i>, quantum phase transitions may occur when K 
changes, and vice versa. 



III. SQUARE LATTICE 

In Fig. |21 we plot the TKE per site versus the flux 
parameter <fi for several values of electron filling i/'s. Due 
to the particle-hole symmetry in the square lattice, the 
DOS always has the property D(uj) = D(—qj), and we 
have E kin (<j),v) = Eki n (4>A - v). Hence E kin -4> curves 
are plotted only for v < 1/2. 

At v = 1/2 (Fig. 0X0, SF^ xy5 has lower TKE than 
the Fermi-liquid phase ((f) = 0) and all other extended- 
SF phases for a given (/>. That is why SF^ x ^ can be 
stabilized befor o 1 ! 16 near half filling. While at v = 3/8 
(Fig. |2Jb)) and even smaller filling factors, SF^ x ^ has 
higher TKE than all other extended-SF phases for a given 
0. At v = 3/8, for < < 0.83, 0.83 < < 1.51, and 
1.51 < < 3.14, SFS xvV SF^ x2V5 and SF? x4 have 
the lowest TKE, respectively. 

The nesting property of the Fermi surface in the 
Fermi- liquid phase at v = 1/2 makes the system resem- 
ble a one-dimensional Perierls system in which a lat- 
tice distortion lowers the TKE of electrons by open- 
ing a gap near the Fermi energy. Therefore the ad- 
vantage of extended-SF phases versus the Fermi-liquid 
phase in TKE can be considered as the generalized 
Perierls orbital-antiferromagnetic instabilities^ of the 
Fermi-liquid phase near the nesting filling, similar to 
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FIG. 2: (color online). Square lattice: total kinetic energy 
per site (in units of t) of extended-SF phases versus the flux 
parameter <j) (in units of <j>o/2n) for several values of electron 
filling i/'s. (a) v = 1/2; (b) v = 3/8; (c) v = 1/3; (d) i/ = 1/4. 

the instability of lattice electrons in a uniform magnetic 
fields. While away from the nesting filling, the advan- 
tage of extended-SF phases in TKE is weakened gradu- 
ally. As for v < 1/4, nearly all extended-SF phases have 
higher TKE than the Fermi-liquid phase. 



and ^2^/2x2V2 com P e ^ e an d show several tricritical 
points. At K > 0.06i, only SF^- x ^ competes with the 
Fermi-liquid phase near half filling. If more complicated 
extended-SF phases are taken into account, the phase 
space will be divided into smaller portions, however, the 
general features should not be altered fundamentally. 

The previous numerical studies£*£*£*2i suggest that the 
SF^ x ^ phase is a key ingredient of the ground state 
in the models of strongly correlated electrons, at least in 
the low doping regime (doping x < 10%, i.e., 0.45 <v< 
0.5), is in agreement with our results. And our results 
also suggest that when v < 0.425 (doping x > 15%), the 
^!/2xi/2 phase does not predominate; when 0.425 < v < 
0.45 (doping 10% < x < 15%), stripe-like extended-SF 
phases, SFf x2 and SF 4x4 should be taken as additional 
competitive ingredients of the ground state. Of course, 
in order to give quantitative evaluation of the TKE or 
the total free energy of interacting electrons, the hopping 
integral t should be renormalized by many-body effects. 

It has been shown before that the instabilities of the 
Fermi- liquid phase in the t-(j> model at nonzero K/t can 
be related to those in the large- N SU(JV) t-J model at 
nonzero t/J 17 . Therefore the phase diagram in Fig. 
which shows different Perierls orbital-antiferromagnetic 
instabilities at different electron fillings can be used to 
predict the existence of similar instabilities in the latter 
model. 



IV. TRIANGULAR LATTICE 




FIG. 3: Square lattice: phase diagram in the v-K parameter 
space (K is in units of i). 

The phase diagram in the v-K parameter space is 
shown in Fig. We notice that SFf x4) SF^ x2V ^ 

and SF^ x3 ^ occupy no space in the phase diagram. 
At v < 0.21, the Fermi-liquid phase is robust against 
all extended-SF phases considered. When v increases, 
SF 4x4 appears in the K < 0.03i region. For 0.3 < 
v < 0.425 and K < 0.06<, SFf x2 , SFf x4 , SF^^ 



For the triangular lattice, the DOS has the symmetry 
D(u,w - (p) = D(-uj,(f>) 19 , consequently, E kln {4>,v) = 
Ekin{^ — <j>, 1 — v). Around the nesting filling v = 3/4, 
the extended-SF phases also exhibit the generic advan- 
tage in TKE versus the Fermi-liquid phase in broad v-cj) 
parameter space (see Fig. Q). 

To obtain the phase diagram as shown in Fig. |5j one 
should notice that the "magnetic" energy per site is K cf> 2 
in the triangular lattice (in which the number of plaque- 
ttes doubles the number of sites). SF^^^ does not ap- 
pear in the phase diagram. For v > 0.82 and K < 0.26t, 
SF^^ dominates the phase space. For K < 0.08i and 
v close to the nesting filling v = 3/4, SFj^, SF^ - 

and SF^ x3 compete. When v is far away from the nest- 
ing filling {y < 0.27), the Fermi- liquid phase is robust 
against all extended-SF phases considered. For larger K 
(K > 0.26t), only SF^^^ is still robust against Fermi- 
liquid around v — 3/4. It suggests that in parallel with 
the importance of SF^ x ^~ near half-filling in the square 

lattice, SF^ x ^ ^ should be a key ingredient of the ground 

state in some lattice-electron models^ in the triangular 
lattice near v = 3/4. 
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SUMMARY AND CONCLUSIONS 



0.0 0.5 1.0 1.5 2.0 2.5 3.0 



0.0 0.5 1.0 1.5 2.0 2.5 3.0 



FIG. 4: (color online) . Triangular lattice: total kinetic energy 
per site (in units of t) of extended-SF phases versus the flux 
parameter <f> (in units of (j>o/2ir) for various values of electron 
filling i/'s. (a) v = 1/3; (b) v = 3/8; (c) v = 1/2; (d) v = 5/8; 
(e) v = 3/4; (f) v = 7/8. 
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FIG. 5: Triangular lattice: phase diagram in the v-K param- 
eter space (K is in units of i). 



In numerous works, the SF^ x ^ phase competes 
with other better-known phases, such as the phases of 
charge-density waves, spin-density waves, superconduc- 
tivity, stripes and Fermi liquid in the square lattice. The 
possibilities of other extended-SF phases in both the 
square lattice and the triangular one, and the compe- 
tition among these extended-SF phases and the Fermi- 
liquid phase, which have never been considered before, 
are the main contributions of this report. 



In conclusion, there are several points to be addressed: 
(a) the generic and strong electron-filling-dependent sta- 
bilization of extended-SF phases in some parameter space 
of the simple t-<p model is an illustration of general- 
ized Perierls orbital-antiferromagnetic instabilities of the 
Fermi-liquid phase in 2D lattices; (b) in the square lat- 
tice, the SF^ x ^ phase is dominantly robust close to 

half filling, which agrees with previous studies; (c) away 
from half filling, other extended-SF phases are possible 
as a result of advantage in TKE against the Fermi-liquid 
phase; (d) in the triangular lattice, we find the SF^ x ^ ^ 

phase is robust around the nesting filling v = 3/4 (in 
a manner similar to the SF^ x ^ phase in the square 

lattice), and other extended-SF phases are also possi- 
ble when v is away from 3/4; (e) since lowering TKE of 
electrons plays an decisive role in minimizing the total 
mean-field free energy in many cases, we keep cautious 
but optimistic hopes about applying our nontrivial new 
results based on a simple model to other more realistic 
systems such as the t- J model and its large- TV limit, the 
Hubbard or extended Hubbard model and their weak- 
coupling limits or large-TV limits, the double-exchange 
model and models with ring-exchange interactions. 
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